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Abstract 

The low-energy amplitude of Compton scattering on the bound state of two 
charged particles of arbitrary masses, charges and spins is calculated. A case in 
which the bound state exists due to electromagnetic interaction (QED) is considered. 
The term, proportional to u> 2 , is obtained taking into account the first relativistic 
correction. It is shown that the complete result for this correction differs essentially 
from the commonly used term Aa, proportional to the r.m.s. charge radius of the 
system. We propose that the same situation can take place in the more complicated 
case of hadrons. 
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1 Introduction 



The electromagnetic polarizabilities a and (3 are fundamental characteristics of the bound 
system. Their magnitudes depend not only on the quantum numbers of the constituents, 
but also on the properties of the interaction between these constituents. Therefore, the 
experimental and theoretical investigation of the electromagnetic polarizabilities are of a 
great importance. In particular, their prediction and the comparison with experimental 
data may serve as a sensitive tool for tests of hadron models. Correspondingly, a large 
number of researchers have been attracted by this fascinating possibility. The electromag- 
netic polarizabilities can be obtained from the low-energy Compton scattering amplitude. 
In the lab frame the amplitude of Compton scattering on the compound system of total 
angular momentum S = 0, 1/2 up to 0(uj 2 ) terms reads |1|, [2|] 

T = T Born + a^i^ei ■ e\ + /3(ki x e x ) ■ (k 2 x e*) , (1) 

where cjj, k, and are the energy, momentum, and polarization vector of incoming (i = 1) 
and outgoing (i = 2) photons (h — c — 1). The contribution T Born corresponds to the am- 
plitude of Compton scattering off a point-like particle with spin, mass, charge, and magnetic 
moment equal to those of the compound system. For spin S > 1 the 0(uj 2 ) part of the 
Compton scattering amplitude contains additional terms, proportional to quadrupole and 
higher multipoles of the bound system ||. In particular, for S = 1 there is a contribution, 
proportional to the quadrupole moment operator. 

The investigation of electromagnetic polarizabilities is interesting not only for systems, 
bound by the electromagnetic interaction, like atoms, but also for those, bound by strong 
interaction, such as atomic nuclei Q or hadrons |J. At present there are many different 
approaches used for the description of the electromagnetic polarizabilities of hadrons: the 
MIT bag model || 0], the nonrelativistic quark model ||, [U], [II], [12[], the chiral quark 
model |14]1 ) the chiral soliton model [15|, [HI and the Skyrme model [[17], [18| . Here 



we mentioned only a small part of the publications on these topics (see also review ||19|| ). 
Though much effort has been devoted to these calculations, all of them can not be consid- 
ered as completely satisfactory. In particular, there is a problem in the explanation of the 
magnitudes of proton and neutron electric polarizabilities within a nonrelativistic quark 
model. It was derived many years ago 0, |20, ^IJ that a can be represented as a sum 

2 _ |HD|0)| 2 . 

where D is the internal electric dipole operator, |0) and \n) are the ground and excited 
states in terms of internal coordinates, and E n and E the corresponding energies. The 
term Aa in a has a relativistic nature contained and its leading term is equal to 

Aa = ^| , (3) 
3M' v ' 

where e and M are the particle charge and mass, is the electric radius defined through 
the Sachs form factor G^. The calculation of the quantity a Q in the nonrelativistic quark 

1 For convenience, our definition of r\ absorbs a total charge e of the system. 
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model without relativistic corrections taken into account leads to the same magnitude of 
a D for proton and neutron. Since Aa is equal to zero for the neutron but gives a significant 
contribution to a for the proton, one has a contradiction between the theoretical prediction 
of a for nucleons and their experimental values, since the latter are close to each other. In 
fact, this approach is not consistent, because there are relativistic corrections to a which 
are of the same order as Aa. 

Starting from second-order perturbation theory one gets the following expression for 

a D : 

a ° 3 -£„)=' W 

where J is the internal electromagnetic current. Using the identity J = i[H, D], where H is 
the Hamiltonian, one comes to the form Eq. (EJ) for a . The relativistic corrections to a in 
Eq. (Q) come from the corrections to wave functions and energies of the ground and excited 
states, and correction to current J. Note, that in ]22| (see also fl2"3| ) it was mentioned, that 
such corrections could exist, but no explicit calculations were made for a realistic system 
and the importance of these corrections was not realized. 

Due to the relation between J and D it is clear that there also is a relativistic correction 
to the electric dipole moment operator (see below) which is connected with the appropri- 
ate relativistic definition of the center-of-mass coordinate. The neglect of this relativistic 
correction leads to an incomplete expression for a, and the missing piece which is calcu- 
lated in the following turns out to be very essential. We expect, that the inclusion of all 
relativistic corrections allow one to remove the big difference between the predictions of 
the nonrelativistic quark model for proton and neutron electric polarizabilities due to the 
difference in Aa. 

The expression for the magnetic polarizability in the nonrelativistic quark model with 
no exchange and momentum-dependent forces has the form PTl E3L EH]. 



- 2 ](n|M]Q)l 2 ef(rf) (D 2 ) 

1 3 ^ E n -E ^ 6m, 2M ' 1 ' 

where M is the internal magnetic dipole operator and the summation in the second term 
on the r.h.s. is performed over the constituent quarks, being the corresponding internal 
radius vector. 

In order to understand the importance of the different relativistic corrections for the 
polarizabilities, it is useful to consider the example of a system where the relativistic cor- 
rections can be obtained analytically. In this paper we calculate the low-energy Compton 
scattering amplitude for a system of two particles with masses mi )2 and charges ei^, bound 
by electromagnetic forces. We consider the case C 1 which provides the validity of the 
nonrelativistic expansion. We consider in detail the cases of spin and 1/2 of the particles 
and give the result for general case of arbitrary spins. 
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2 Scattering amplitudes 



For the electromagnetic interaction between particles a simple estimate shows that the part 
t of Compton scattering amplitude, proportional to uj 2 has the form: 

t = u>^L+c 2 ^ + cA 2 +--\ , (6) 

where a = l/((ig) is the Bohr radius, g = — eie 2 > 0, \i = m 1 m 2 /(m 1 + m 2 ) is the reduced 
mass, Eq = —fig 2 /2 is the ground state binding energy in the nonrelativistic approximation, 
and q are some quantities, bilinear with respect to e 1; e 2 and depending on the ratio of 
charges and masses. Here, for defmiteness, we assume u — u±. The two first terms of this 
expansion contain the parameter g <C 1 in the denominators and, therefore, come from 
the contribution of big distances r ~ a (or small momenta p ~ 1/a — fig) to the matrix 
element. These two terms which have no contributions from the Born amplitude are the 
ones we are going to calculate in this article. Since they are determined by a contribution 
from big distances (small momenta), it is possible to use the nonrelativistic expansion in the 
calculations. In fact, the first term is known and is proportional to the electric polarizability 
a Q , calculated in the leading nonrelativistic approximation (see below). Some contributions 
to the second term are also known, namely, those, containing the magnetic polarizability 
(3 Eq. (|5|) and the correction Aa Eq. (|3|). These contributions come from the expansion 
of the photon wave functions over kr ~ uj/fig and from the expansion of the propagator 
of the system with respect to the photon energy and the center-of-mass kinetic energy of 
intermediate states. The corresponding results for j3 and Aa can be obtained using the 
nonrelativistic Hamiltonian of the system. As was mentioned above, the other source of 
the contributions to c 2 , which has not been investigated so far, is the relativistic correction 
to the Hamiltonian of the system and the corresponding corrections to the wave functions, 
energy levels, and currents. We will obtain the complete result for the second term in 
the expansion Eq. (^). In the expression for the electromagnetic current we neglect for a 
while the dependence of the form factors on the momentum transfer. We will take this 
dependence into account at the consideration of the general case of arbitrary spins. 



2.1 The system of two spin-0 particles 

Let us consider first the bound state of two spin-0 particles. In order to calculate the 
Compton scattering amplitude, it is convenient to put the system into the external elec- 
tromagnetic field A(x, t). In this case the nonrelativistic Hamiltonian has the form 

H nr [A] = -gL + -gt - - 9 , , (7) 
2mi 2m 2 |ri — r 2 | 

where 7Tj = pj — ejA(rj,t). Let us pass to the variables r and R, corresponding to the 
relative and center-of-mass coordinate: 

m 2 mi 
r x = R+— r, r 2 = R-— r, M = m 1 + m 2 . (8) 



Then, the momenta are 



p^p + p, p^^p-p. w 
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where P 



-zVfl and p - 

Hri.r \ A 



iV r . For A = we have 

p2 



0] 



2M 



+ nr 2M + 2/i r " 



(10) 



The first relativistic correction Hb (Breit Hamiltonian, see, e.g., f2"6fl ) to Eq. (^) reads 



H R \A] 



2^ 2 



+ 



7T 



m 3 . 



+ 



9 



2mim 2 



— + — ? 



2 



7Ti 7T 



l"2 



The first term in Eq. (|TT|) is the correction to the kinetic energy and the second one is the 
correction due to magnetic quanta exchange, corresponding to the space component of the 
photon propagator in the Coulomb gauge. If A = then in the center of mass frame where 
the eigenvalue of the operator P is equal to zero we have 



#b[A = 0] 



p=o 



H 



B 



8 V m? mo / v 
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2m\m2 



(12) 



The terms, containing the operator P in the Hamiltonian, determine the contribution of 
recoil effect to the Compton scattering amplitude. Within the precision of the present 
calculations these terms should be taken into account only in the Hamiltonian H nr and can 
be ommited in Hb (see below). The correction 5eq to the ground state energy, related to 
the Hamiltonian Hb reads 



Se =(0\H B \0) = -g 4 



5 4 / 1 n /i 3 

8 ymf m^l mim 2 



(13) 



Let us start the calculation of the Compton scattering amplitude with the amplitude T nr 
obtained with the use of the nonrelativistic Hamiltonian Eq. fllPl) . This amplitude can 
be represented as a sum T nr = T res + T s of resonance and seagull parts. The part T res is 
determined by the second order of perturbation theory with respect to the terms in H nr [A], 
linear in the vector potential A. In the lab frame it has the form 



T 



■(ip \ exp[-z(ki 

p2 



k 2 )Rle* • 



— piexp(-ik 2 ri) + — p 2 exp(-ik 2 r 2 ) 

mi m 2 



x 



x [e + uj 1 - 

+ (6! 



2M 
e* 2 , ut <- 



— H nr ] 1 e 1 ■ 



-w 2 , ki 



— piexp(zk 1 r 1 ) + — p 2 exp(zkir 2 ) 

m\ m 2 



|V>o> + 



(14) 



Here ipo(r) — n^ 1 ^ 2 (ng) 3 ^ 2 exp(-figr) is the wave function of the ground state, depending 
on the relative coordinate r. The final momentum of the bound system is equal to k x — k 2 . 
Using the relations Eq. (|8[) and Eq. (|9|) and making a simple transformation in order to 
cancel the exponents containing R, we obtain 



T 



•(^o|e 2 



ei mi m 2 e 2 m 2 m x 

-(P + exp(-,-k 2 r) - _(p - _ kl ) exp(,-k 2 r) 



x 



x G(^i) ei • p 
+ (ei <-> e* 2 , cji 



ei ,.m 2 e 2 .mi 

— exp(«— k x r) exp(-z— k ir ) 

mi M m 2 M 



-tx>2 , ki <-> — k 2 ) 



(15) 
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Here G{uj) = [e + u> — u 2 /2M — H nr ]~ l is the nonrelativistic propagator of the system in 
the operator form. The seagull amplitude T s is determined by first order of perturbation 
theory with respect to the terms in if nr [A] which are quadratic in A. Similar to Eq. (HI 
we obtain 



T s = -e 1 ■ e^Vol 



— exp(^(ki - k 2 )r) + ^_ eX p(-z^(k 1 
rrii M m 2 M 



.rrii 



k 2 )r) 



Wo) ■ (16) 



Performing the expansion of Eq. (|14D and Eq. ([TBD with respect to k 12 and cj 12 up to 
quadratic terms and using the relation uji — u 2 = (ki — k 2 ) 2 /2M , we obtain 



T 



" ei€2 M + 
-[ei x ki][ea x k 2 



eie 2 u) 
2?' 



2/V 



gi _ _e2_\ 2 gi + e 2 / < , 
mi m 2 / Mg 2 



+ 



mv 



m 2 



+ 



2Mg< 



ei 
mi 




e 2 
mi 



(17) 



There is no need here to distinguish between uii and u; 2 in the 0(oj 2 ) term. Therefore, we 
set u)i = oj 2 = u) in Eq. (fl7|) . The result Eq. (|17D is in agreement with Eq. (§) and Eq. @, 
with a D calculated in the nonrelativistic approximation, since in our model 



ei 



e 2 



mi m 2 



«onr = o(V'o|D nr <3oD fM .|V'o) 
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3 

ei + e 2 
3M 



gi e 2 y 
2/zg 4 Vmi m 2 / 



(V'oleir 2 + e 2 r||^ ) 



ei + e 2 / ei e 2 
Mg 2 Imf m 2 



<lfo|M|Vwo> = , (^o|D 2 |^o) = ^ f — - — ) 



g 2 Vmi m 2 / ' 



(V>o 



2^,2 



err 



1' 1 



+ 



ejr 2 



g 2 Vmf 



+ 



m 2 



mi m 2 

Here Go is the reduced Green function in the operator form: 

Go = [eo-H nr + tO]- 1 (l-\^ }(^ \). 



(19) 



The details of calculations of different matrix elements, containing the operator Go are 
presented in Appendix A. 

We pass now to the calculation of the relativistic corrections to the electromagnetic 
polarizabilities, connected with the Breit Hamiltonian Hb Eq. ([0]). We perform the cal- 
culations in the same way as at the derivation of Eq. fll4|) and Eq. (]T6|), but for the Hamil- 
tonian H = H nr + Hb- At the calculation of these corrections within our accuracy the 
terms of the Hamiltonian Hb[A], quadratic in A do not contribute to the electromagnetic 
polarizabilities, i.e. the seagull contribution from Hb[A] is absent. In the corrections to 
the electromagnetic polarizabilities from the resonance part of the amplitude, we can take 
the second order of expansion with respect to u> of the operator Green function and put 
ki )2 = elsewhere. This means that within our accuracy we can neglect in Hb[A] the terms 
containing the total momentum P and replace the exponents in the photon wave function 
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by unity. Therefore, the terms in the Hamiltonian Hg[A], linear in A, can be represented 
in the form — A(0)J#, where 



g( e i - e 2) ( + JLr r . \ | ( 9 0) 

m\ irio) 2 2mim 2 r V r 2 ' 



is the correction to the operator of the total internal current J : 

J = J nr + J B = (ei/mi - e 2 /m 2 )p + J B • (21) 



Let us now discuss the relativistic correction to the electric dipole moment operator. In 
the lab frame it is equal to ~D tot = e^i + e 2 r 2 . For the total Hamiltonian of the system the 
following relation holds 

H„ [A = 0] + H B [A — 0], D w ] = (22) 
= ^(l-^)-^ p 1 + 4(r. Pl ))+(l-2). 

The r.h.s. of Eq. (|22[) is nothing but the operator of total current in the lab frame, con- 
taining the total momentum P. This can be verified by differentiating the Hamiltonian 
H nr [A] + Hb[A] over A at A = 0. Therefore, there are no relativistic corrections to the 
total electric dipole operator D tot . The internal electric dipole moment operator is defined 
as 

D = D toi - (ei + e 2 )R cm , (23) 

where R cm is the center-of-mass vector. This vector is defind in such a way that it satisfies 
the following relations: 

[R cm , P] = % , i[H tot , R cm ] = — - , (24) 

n-tat 

where H to t is the total relativistic Hamiltonian of the system, and P is the total momentum. 
Within our accuracy the second relation in Eq. (^4j) reads: 

i[H nr \A = 0] + H B [A = 0] , R cm ] = ^ (l - Hnr[ M =0] ) ' (25) 



It is known (see, e.g. , P7| ) that there is a relativistic correction to R cm in classical 
electrodynamics. For the case of two particles, interacting due to electromagnetic field, the 
corresponding operator which satisfies the relations Eq. (123) has the form 



R cm = R+— (L,tf_-°-\ + L, .pL- a\_( r , vL + vL_d) . (26) 

2M\\ 1, 2m l 2rJ \ ' 2m 2 2rJ \ ' 2m 1 2m 2 rjj V ; 

Here we took into account the first relativistic correction and use the notation {a, b} = 
ab + ba. In terms of the variables r and p (see Eq. (Q) and Eq. (|9])) we obtain : 

R cm = R + (m2 2M2 m0 ({r, H nr } + g r ~\ , (27) 
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where the term proportional to the total momentum P is omitted. Substituting this ex- 
pression into Eq. (|23| ) we obtain the relativistic correction to the internal electric dipole 
moment: 

(ei + e 2 )(mi - m 2 ) (, „ , r 



D 



B 



2M 2 



{r, H nr } + g- 
r 



Note that, as should be the case, the operator of total internal current J = J nr 
within our accuracy the relation 



i[H nr + H B , D 



D 



B\ 



(28) 
+ Jb satisfies 

(29) 



where D nr is defined in Eq. (pT8|) , and H nr + Hb is the internal part of Hamiltonian (see 
Eq. (|10D and Eq. (0)). 

Let us write down now the corrections to the 0(u 2 ) term of Compton scattering ampli- 
tude, related to the Breit Hamiltonian and the corresponding current. The correction due 
to J B reads 



(30) 



The 0(u 2 ) correction to the amplitude, connected with the expansion of the propagator 
with respect to Hb, has the form 



(31) 



-io 2 (^o\e* 2 ■ J nr [GqHbGq + G s H B Go + G H B G$e x ■ J nT # > + (ei <-> e^) . 
The contribution due to the correction to wave function is 

t w = -u 2 (ipo\[^*2 ■ JnrGlei ■ J nr G H B + H B G e* 2 ■ J^-G^ei ■ J nr ]\i/}o) + (ei <-> e* 2 ) . (32) 
At last, the contribution corresponding to the correction to the ground state energy reads: 



t e = 3u 2 6e {i/>o\eZ ■ J„ r Go £ i ' J n»#o) + (ei <-> e 2 



In order to calculate the matrix elements in Eq. 
following relations (see Appendix A): 



(33) 



-Eq. (|33[ ) it is convenient to use the 



G rip = 




r 

2" 










( r 3 


G r 2 ip = 


















G r 3 ^ = 











ar - 

5ar 3 



11a 3 



2 

5a 2 r 2 



r 

9 



155a 4 



r 5ar 5a \ , 

3 + ^ + -r c 



(34) 



^0 



where a = l///a. Using (|34]) and also the relation p = ifi[H nr ,r 
results for the contributions Eq. (|30|)-Eq. (|3 



2/j,u; 2 e i 



9g 2 

H 2 (jj 2 ei ■ e 2 



H 2 uj 2 Ei ■ e* 2 



H 2 uj 2 ei ■ e* 2 



ei 
mi 



e 2 
m 2 

gi _ e 2 y 

mi m 2 J 
ei _ e 2 y 
mi m 2 / 

ei _ _e2_ 
mi m 2 



we get the following 

89 (ei-e 2 )~ 



m 



rrin 



1061 / 1 



288 



4 I mf 



^mj 

1 1 

+ — 

m 



m| 



+ 



+ 



4 yumim 2 
25 
3/imim 2 
14 



(35) 



2 129 



1 1 



mv 



mf 



9yumim 2 
+ 



yumim 2 



Representing the sum of all contributions in Eq. 
result for a oB : 



5D as ijj €i ■ e\(x B we obtain the following 



OLoB = —- 



1 

7 




113 \ (ei + ea) 



4M 



mi 



e 2 
m 2 



[mi - m 2 J 
2mim 2 



(36) 



If one starts the calculation of ct Q from Eq. @, then it is easy to check that the first 
term in Eq. fl3"E|) corresponds to the sum of corrections due to modification of wave func- 
tion, propagator, and ground state energy. The second term in Eq. (36) corresponds to 
the contribution due to the relativistic correction to the electric dipole moment Eq. (p8|). 
Therefore, this correction appears due to the correct description of the center-of-mass mo- 
tion. It is seen that the first term in Eq. ([36]) has the same dependence on charges as a onr 



while the second term is similar to Act , Eq. ([Tj|). Taking a sum of ct B , ctonr, and Act, we 
come to the following result for ct for the system under consideration: 



a 




(37) 



Thus, the relativistic corrections to ct has reduced to a renormalization of ct onr and essen- 
tially to a modification of Act , Eq. ([3D. One can expect that the last statement is valid not 
only for the system under consideration. Indeed, due to the definition of D, the correction 
to at related to the modification of R cm is proportional to Q/M , where Q is the total 
charge of the system, and, therefore, has the same structure as Aq. 

As a nontrivial test of our method of calculation we checked the fulfilment of the low- 
energy theorem for the Compton scattering amplitude. At u = this amplitude should 
have the form 



T(u = 0) 



(gi + e 2 ) 

En 



-61 



(ei + e 2 y 
M + e 



-ei 



(ei + e 2 ) 
M 



£o 
M 



(3i 



where Eq = M + En is the mass of the system. It is interesting, that the term in r.h.s. 



of Eq. (]380 , proportional to the nonrelativistic energy e 0) appears as a contribution of 
terms from the Breit Hamiltonian Hb[A], which we checked by explicit calculations (see 
Appendix B). 



2.2 The system of a spin-0 particle and a spin-1/2 particle 

Let the first particle have the spin 1/2 and the second particle have the spin 0. Then we 
should add the term 

5H nr [A] = - 6l(1 + /tl) si • H (39) 
mi 

to the nonrelativistic Hamiltonian H nr [A], Eq. ([/]). Here H is the external magnetic field, 
Si = <Ti/2 is the spin operator of the first particle, and K\ is its anomalous magnetic 
moment in units ei/2m-i. There is also some additional contribution 5Hb[A] to Hb[A], 
Eq. ([12]) (see, e.g. f26|). The terms of 5Hb[A] linear in si as well as Eq. ( |39"]) determine the 



S 



O(oj) terms of the Compton amplitude. These terms are well-known and follow, together 



with the a;- independent term, from the low-energy theorem [ 2!| . As it was explained in the 
previous subsection, it is sufficient within our accuracy to account for the Breit Hamiltonian 
only in the long-wave limit, i.e. at uo\^ = in order to obtain the 0(lv 2 ) terms of Compton 
amplitude. In this limit the Hamiltonian 5Hb[A] reads 

5H B [A] = - 21 2 l> tt5 r + "T s i • (r x tti) ) + — ^ • r x tt 2 . 40 

2mf V r 6 / mim 2 r 

The explicit calculation shows that the contribution of SH nr [A] given by Eq. (^) and the 
terms in Eq. ( |40|) linear in o"i do not lead to any contributions to c lj2 in Eq. (|[), i.e. they 
can be neglected in the calculation of polarizabilities within our accuracy. In particular, 
there are no terms linear in the spin in the quantities Ci j2 , which is in agreement with the 
general conclusion on the absence of terms 0{uj 2 ) linear in spin in the non-Born part of the 
Compton amplitude p2 |. The only term which should be taken into account in addition 



to those considered in the previous subsection, is the spin-independent term in Eq. ([JO) 
(Darwin term): 

^Ml±^ (r , (41) 

It follows from Eq. (|40f) that there is no correction to the current associated with the 
Hamiltonian 5dHb- Using the expressions ( |31| - |3"3|) with the replacement Hb — > SdHb and 
the relations Eq. (|4j) we obtain 

S D t e = 0, 6 D t w = - ^fl 1 f 2 (il - ^-Y (1 + 2k,) , (42) 

d D t p = 0, d D t e = „ (1 + 2%). 



im 2 g 2 Vmi m 2 / 

As a result, the correction to the electric polarizability associated with the Breit Hamilto- 
nian in the system of spin-0 and spin- 1/2 will be the sum of a B Eq. ( |5o| ) and 

*a ofl = ^ - -^V (1 + 2*) . (43) 
2mfg z \mi m 2 / 



2.3 The system of two spin- 1/2 particles 



In the case of two spin- 1/2 particles it is necessary to account for two Darvin terms in 
addition to the Breit Hamiltonian Eq. (O), corresponding to both particles 



7Tfl(l + 2Ki) 7rg(l + 2/t 2 ) 

^ = 2m? 5(r) + 2ml ^ 

and the Hamiltonian, corresponding to spin-spin interaction 



(44) 



5M 



0(l + «i)(l + K2) 



B 



mim 2 



3(n- Si)(n- s 2 ) - Si ■ s 2 87r 
= h —0(t) Si • s 2 



(45) 
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with n = r/r. It is more convenient to rewrite 5 s Hb in terms of the total spin operator 
S = Si + s 2 : 



g(l + K 1 )(l + K 2 ) 
o s H B = - 



SniTijQij . . / 2 



2r 3 



+ 4n5(r) (-S 2 - l) 



(46) 



2mi7B2 

where the operator Q^, quadratic in S, is equal to 

2 

Qij = SiSj + SjSi - -S 2 5ij . (47) 

Note that in such a system as positronium it is necessary to add the contribution of the 
annihilation diagram, which results in the replacement (2S 2 /3 — 1) — > (7S 2 /6 — 1) in the 
coefficient of the ^-function in Eq. (flff ) (of course, in this case mi = m 2 , e\ = — e 2 ). As in 
the previous subsection, the terms proportional to the ^-function in Eq. ( f4"4"D and Eq. (fK)|) 
give the contribution 8a Q B, which should be added to a OB , Eq. (j36l): 



31/i ( e 1 _ e 2 V 



2g z \rrii m 2 J 



+ 1 + + 4(1 + K-QU + K2 ) ( 2 S{S + l} _ l 



mi m 2 mim 2 \3 

(48) 

Here we replaced S 2 by its eigenvalue S(S + 1), where S = 0, 1 is the total spin of the 
system. The term in Eq. (48) containing the tensor operator 3(n ■ S) 2 — S 2 determines the 



contribution to the 0(uj 2 ) part of the Compton amplitude, which has the form 

t(tensor) = uj 2 a T e\e 3 2 * (Qij) , (49) 

where (• • •) denotes the averaging over the spin part of the wave function. Of course, 
t(tensor) vanishes if S = 0. Since there is no correction to the current or to the energy of 
the ground state due to the tensor part of 5 s Hb, the contributions to ax come only from 
the corrections to the propagator and to the wave function. Using Eq. Q3"1"D, Eq. (|3"2]), and 
the relations (see Appendix A) 

Go^.-r 2 ^) = - ^-J Hii + |Vb) (50) 
G (3r irj - r%j)r\^) = - ^ ~ ^ ( ^ + ^ + ^ [^o) , 

we obtain 

47(1 + «i)(l + k 2 ) ( ei e 2 \ 2 
" T 40Mc/ 2 Ui m 2 y 1 J 

In the system of two spin-1/2 particles there is a big paramagnetic contribution to the 
magnetic polarizability from the first term in Eq. (|5|). The main contribution corresponds 
to the transition from the ground state with the total spin S — to the state with S — 1, 
with both states having the same angular momenta, / = 0, and radial quantum nambers 
n r = (hyperfine splitting). Representing the spin part of the magnetic moment operator 
in the form 

TV/T £ £ £ e «(l + K i) 

M s = fisi + f 2 s 2 , fi = , 

rrii 
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and using Eq. fl45|) , we obtain: 

A = -%^ , -=^. (52) 

As was pointed out in the previous section, for positronium it is necessary to change the 
coefficient of 5-function in Eq. ( f46|) . As a result, the contribution of the first term in Eq. (|5|) 
to the magnetic polarizability of positronium is: 

Pi = ±^a 3 , (53) 

where upper sign corresponds to parapositronium (S = 0), and lower sign to orthopositro- 
nium (S = 1). 



2.4 The system of two particles with arbitrary spins 

Let the particles have the spins Si^ and magnetic moments which we represent in the 
form 

/x a = ^(l + K a ), a=l,2. (54) 



The electromagnetic current for each particle has the form (see, e.g., p9| , |30[| ) 



2m 2m 



^(p) , (55) 



where q = p' — p. The operator E Ml , is a generalisation of the corresponding matrix for spin 
1/2. The indices numerating the particles have been omitted. The quantities F e and G m 
depend on q 2 and (s^q 11 ) 2 , where s M is the 4-vector of the spin operator. This quantities 
are normalized as follows: 

F e (g = 0) = l, G m (q = 0) = l + (56) 

If we neglect the g-dependence of the form factors, then, in addition to the Breit Hamilto- 
nian for two spin-0 particles, Eq. fllT]), it is necessary to take into account the Hamiltonian 
Eq. (J45[) (with the corresponding spin operators) and two other contributions | 29fl . Namely, 
the Darwin Hamiltonian 

5 D H B = £ + ^a){s a + Ca)S(r) , (57) 

0=1,2 dm a 

C = for integer spin and £ = 1/4 otherwise, and the term containing the quadrupole 
moments of the particles: 

SqHb = E " (1 + f 3 a) S (n • s a ) 2 - s 2 a ) , (58) 

£ = 1/ (2s — 1) for integer spin and £ = 1/ (2s) otherwise. It is clear that all matrix elements 
can be calculated in the same way as in the previous subsection. The averaging over the 
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spin variables can be done using the following relations 



(SS Z \ 
(SS, 



SliS 2 j + S 2 iS lj - -(%SiS 2 



SliSlj + SuSij SijSi 



2 
3' 



S,S' Z ) = A(S,s 1 ,s 2 )(SS z \Q ij \S,S' z ), 
S,S' Z ) = B(S,s 1 ,s 2 )(S S z \Qij\S,S' z ) 

~1), and for S > 1 



(59) 



where S = S! + s 2 is the total spin operator, is defined in Eq 

A 2 + 2A(Ax + A 2 ) -3(Ai - \ 2 f 



A(S,s 1 ,s 2 ) 

5(5,51,52) 



2A(4A - 3) 

3A 2 + A(2Ai - 6A 2 - 3) + 3(Ai - A 2 )(Ai - A 2 - 1) 
2A(4A - 3) 



(60) 



Here A = S(S + 1), Ai )2 = 51,2(^1,2 + 1) are the eigenvalues of the operators S 2 and s 2 2 , 
respectively. For S = 0, 1/2 we put A = B = 0. As a result, we obtain the following 
generalization of Eq. (|48D to the case of arbitrary spins: 



5a. 



oB 



62/i / ei e 2 \ 2 
3g 2 \mi m 2 J 

2(l + «i)(l + K 2 



+ 



mim 2 

The generalization of Eq. ( |5"ID is 

47// / ej_ _ ^2_\ 2 
40(/ 2 Wi m 2 / 
1 + 2ki 



1 + 2«i , . 1 + 2^ 

— (Sl + Cl + — («2 + C2J + 

mf m 2 
■ (A — Ai — A 2 ) 



(61) 



2(i + K 1 )(i + « a ; 



A(S,s 1 ,s 2 )+ 



+ 



£iB(S, st, 8 2 ) + 1 + 2 2 K2 &B(S, s 2 , st) 
mi mi 



(62) 



Let us now take into account the g-dependence of the electromagnetic form factors of the 
constituents defined in Eq. (|55|) . We assume, that the scale of variation of these form 
factors are much larger than the typical momentum transfer ~ fig. In other words, the 
characteristic size of each constituent is much smaller than the size of the whole system 
a = 1/fig. In this case it is sufficient, within our accuracy, to take G m = 1 + k and to 
expand the form factor F e up to quadratic in q terms: 



i^(g 2 ,(^) 2 )^l-^! + r " (S ' q)2 



(63) 

6 2 K ' 

where r 2 es are some constants. Multiplying the 0(q 2 ) terms in this expression by — Airg/q 2 
and performing the Fourier transform, we obtain the additional terms in the Hamiltonian 



S f H B 



E 

a=l,2 



2ng 



2 22 

ea sa a 



) S(t) + 



2 3(s a • n) 2 - s- 



2r 3 



(64) 



Since the terms in this Hamiltonian have the same structure as above, it is easy to write 
down the result for the corresponding corrections to polarizabilities: 

62/z ( e\ e 2 x 2 



5fa oB 



3g 2 Vtoi m 2 
47/x / e\ 



(•■2 



40(/ 2 \mi m 2 



r ll ~ rll^l + r e2 - r s2^2) 
2 . 

r 2 sl B(S, s u s 2 ) + r 2 s2 B(S, s 2 , s t ) 
12 



(65) 



If the parameters of the form factors r 2 s ~ 1/m 2 <C a 2 , then the contributions Eq. ( |65f) 
to the polarizabilities are of the same order as a B- The first relativistic correction to the 
Compton scattering amplitude at uj = 0, Eq. (|38|), is proportional to e = —fig 2 / 2 and 
is independent of the spins of the constituents. Then, the correction to the amplitude at 
u = connected with spin-dependent terms in Breit Hamiltonian as well as the Darwin 
terms (also having the spin origin) should vanish. This statement was checked explicitly 
(see Appendix B). 

Let us consider now the paramagnetic contribution to the magnetic polarizability from 
the first term in Eq. (|5]). Let si > s 2 . Then, the total spin of the ground state is S = si — s 2 , 
and the main contribution corresponds to the transition from the ground state to the state 
with S — Si — S2 + 1, with both states having the same angular momentum, I = 0, and 
radial quantum namber n r = (hyperfine splitting). A simple explicit calculation leads to: 

4/i/2(si - s 2 + l) 2 

This term should be added to the diamagnetic contribution j3di a (see Eq. (|l 

= L ( iL + 3.\ — (fL _ fl ) ( 6 7) 

2g 2 \m\ ml J 2Mg 2 Wi m 2 



3 Conclusion 



We have obtained the complete result for the first relativistic corrections to the electromag- 
netic polarizabilities, including the tensor part which exists for the total spin S > 1. We 
demonstrated that, within our accuracy, this tensor part contains the quadrupole moment 
of the system and no any higher multipoles. For the system of two spinless particles it is 
easy to check that the total relativistic correction Eqs. flTBp and ( |3"6j ) is negative at arbi- 
trary masses and charges. In the general case of non-zero spins and arbitrary anomalous 
magnetic moments the relativistic correction Aa + a B + Sa B, where 5a B is given by 
Eq. (0), can be positive. It is interesting to consider some special cases. The first of them 
is a hydrogen-like ion. In this case e\ = e, e 2 = —Ze, and m 2 3> m\. In the limit m 2 — > 00 
the result for electric polarizabilities is independent of the spin and magnetic moment of 
the nucleus. Neglecting also the anomalous magnetic moment K\ of the electron, we obtain 
from Eq. (0), Eq. (0), Eq. Q3§), and Eq. 
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a = 2m 3 a3 m Z 4 ~ 3m 3 a em Z 2 ' (68) 

where a em = e 2 = 1/137 is the fine-structure constant. Note that in this limit the correction 
Aa , Eq. @, vanishes. The result Eq. ( |68|) is in agreement with that obtained in [31]] with 



the use of the reduced Green function of the Dirac equation for an electron in a Coulomb 
field. For the magnetic polarizability at s 2 = we have 
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For s 2 — 1/2 in the limit m 2 3> m x there is a very big contribution from the paramagnetic 
part of the magnetic polarizability, Eq. fl52|). In the Compton scattering amplitude this 
contribution should be taken into account only for photon energies uj much smaller than the 
energy Ehf ~ «L m i/ m 2 °f the hyperfine splitting. For /ig 2 us Ehf the paramagnetic 
contribution should be omitted. 

Another interesting example is positronium. As we mentioned above, in this case it is 
necessary to replace (2S 2 /3— 1) — > (7S 2 /6 — 1) in the coefficient of the 5-function in Eq. ( |4~6|) 
due to the contribution of the annihilation diagram. Putting m\ = m<i — m, e% — — e 2 = e, 
and Ki = k<i = 0, we obtain Act = and the complete result for the polarizabilities 

36 1 J -1001 forS = 

(ma em ) 3 Qm 3 a em \ 735 for S = 1 

^ ^ ^ 56m 3 a 3 , m m 3 a em ^ 

As in the previous case, for photon energy u 3> OL\ m m the paramagnetic contribution should 
be omitted in the Compton scattering amplitude. 

For 5 = 1 (orthopositronium) we also have the tensor polarizability 

47 

ar = -20^Aw (71) 

Thus, we have shown that the complete set of the first relativistic corrections differs es- 
sentially from the commonly used term Aa. We suppose that for the electromagnetic 
polarizabilities of hadrons investigated within the constituent quark model an analogous 
situation may be found. 
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Appendix A 

In this Appendix we derive the formulas for the result of the action of the operator Go, 
Eq. flT5|), on the wave function \ipo), multiplied by some polynomial of r. More precisely, we 
obtain the expression for GoYi m (r /r)r n \ip ) in the form of the product Yi m (r/r)P(r)\ijjo), where 
P(r) is some polynomial. Since the Hamiltonian H nr commutes with the operator of an- 
gular momentum 1 = r x p, we can make the following transformation: 

G Y lm (r/r)r n \i; Q ) = Y lm (v , W)Gf \r n - ^ (r n ))|^ ) , (72) 
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where = [e - H®. + iO]~\ = -(2/ir)" 1 9 r 2 r + 1(1 + l)/(2^r 2 ) - g/r is the radial 
Hamiltonian with the angular momentum I, and 

{r n ) = {Mr n \^)= K ~ 1 ^a\ (73) 
where a = 1/ [xg. In the derivation of Eq. (|72|) we used the identity 

(1 - |^o)(^o|)^ m (r/r)r n |^o) = Y lm (v/r)(r n - 5 l0 {r n ))fa) . 

For our purposes it is sufficient to consider the cases n > I — 1 for I ^ and n > 1 for 
/ = 0. It is easy to check that in these cases one can represent the result of action of G$ 
in r.h.s. of Eq. ( [72]) in the form 



GS\r n -5 l0 (r n ))\^ ) = ^C k r k fa) 1 (74) 

k=0 



where Ck are some constants to be found. Acting on both sides of this equation with the 
operator Eq — H®. and collecting the coefficients with different powers of r, we obtain 

r--WO = (75) 

+ ^ {k - l+2){k+l : % ^_ g{k+i) , 

k=0 \ A ^ ) 

From this relation we can find the coefficients Cj. For the case n > I — 1, I ^ 0, we finally 
obtain 

PV/MB| , . ^ . . >-Z + l)! (n + Z + 2)! "+ 1 (/c - 1)! {2r / a) k 
G Y lm (r/r)r fa) = -r»»(r/r) g {2/a)n+1 {n+ 1}l E (Jfc + f + 1)l l*>> 

(76) 

For the case n. > 1, I — we have 

(n + 2)! ^l/(2r/a) fc fc + 2\. ; . 

Using the formulas Eq. (f75|) and Eq. ( [76] ) one can easily calculate all matrix elements 
needed. 



Appendix B 

In this Appendix we check the fulfilment of the low energy theorem. Namely, we reproduce 
the two first terms of the expansion with respect to Eq/M of the Compton scattering 
amplitude at to = 0: 

(ei + e 2 ) 2 / e \ 



r(w = 0)»-« 1 .«^-VMl-£). (78) 
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In fact, the first term is contained in Eq. (|17D . In order to obtain the second term, we have 
to take into account the corrections to the current, seagull, wave function, propagator, and 
energy due to the Breit Hamiltonian Hb[A]. The contribution to the amplitude at to = 
due to Jb reads 



■<Vo|[C2 " J bG € 1 ■ J nr + 6* • J nr G € 1 ■ J B ]|^o) + 



(79) 



The contribution due to the correction to seagull (the terms in Hb[A] being quadratic in 
A) reads 

^2" 



<^o| 



+ 



_L _|_ 

m\ m| 



+ 



iO) 



mirri2 



l^o). 



The contribution connected with the expansion of propagator with respect to Hb has the 
form 

T p = -(ip \(el ■ J nr G H B Go e 1 ■ J nj # ) + (ei <-> e 2 ) . (81) 
The contribution due to the correction to wave function is 

T w = -(^o\[e* 2 ■ JnrG e 1 ■ J nr G H B + H B G e* 2 ■ J nr G e 1 ■ J nr \\ipo) + (ej «-> e 2 ) . (82) 

At last, the contribution corresponding to the correction to the ground state energy reads: 

T e = 5e (i)o\e*2 ' Jnr^o 6 ! ' Jnrl^o) + (^1 62)- ( 83 ) 

Using the results of Appendix A, we obtain the following expressions for the corrections: 



T 



T 

-L it 



T 



e l ' e 2 2 2 / e l e 2 

^—9 M 

6 \mi m 2 

5/x I — ~ H : 



5/i 



( e i 



'2 2 

-9 t* 



+ 



mi" 



e i ' e 2 
12 

6i ■ £ 



2 3 / e l e 2 
5 M 

mi m2 



m 2 - 
2 



7/i 



Vmi 3 
8ff 
mi m2 
1 



e 2 
m 2 J 



+ 4 



ei - e 2 



6 



ei • e 2 2 
-5 A* 



ei _ e 2 \ 
mi m 2 / 

3 _ _£2_\ 



9/i 

2 r 



1 



+ 



m 2 " 
1 

m 2 s 



mi m2 



12 



54) 



mi 7712 
14 



4 ' Vmi mi J 
Summing up these contributions, we get 

Tb = — 61 • €o 



5/x 



Vmi 



1 1 

+ 



m 2 - 



mi m,2 
8 

+ 



mi m 2 



[ei + e 2 J 



-MP 



5) 



" 2 2M 2 

which is the second term in Eq. ([78] ) . Let us consider now the contribution to the Compton 
amplitude at uj = 0, connected with the spin-dependent terms and the Darwin terms in 
Breit Hamiltonian. Note that all these terms are proportional to either 5(r) or to the 
operator (3riirij — 5ij)/r 3 . The terms oc 5(r) give the contributions to T w and T e . Using 
the results of Appendix A, it is easy to show that the sum of these two contributions is 
zero. The terms oc (3niTij — Sij)/r 3 give the contributions to T w and T p . Again, direct 
calculations show that they also cancel each other. Therefore, we proved, that the first 
relativistic correction to the Compton amplitude at uj = is spin- independent, which is in 
agreement with the low-energy theorem. 
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